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ABSTRACT 

The excursion set theory, where density perturbations evolve stochasticahy with the 
smoothing scale, provides a method for computing the dark matter halo mass func- 
tion. The computation of the mass function is mapped into the so-called first-passage 
time problem in the presence of a moving barrier. The excursion set theory is also 
a powerful formalism to study other properties of dark matter halos such as halo 
bias, accretion rate, formation time, merging rate and the formation history of halos. 
This is achieved by computing conditional probabilities with non-trivial initial con- 
ditions, and the conditional two-barrier first-crossing rate. In this paper we use the 
path integral formulation of the excursion set theory to calculate analytically these 
conditional probabilities in the presence of a generic moving barrier, including the one 
describing the ellipsoidal collapse, and for both Gaussian and non-Gaussian initial con- 
ditions. While most of our analysis associated with Gaussian initial conditions assumes 
Markovianity (top-hat in momentum space smoothing, rather than generic filters), the 
non-Markovianity of the random walks induced by non-Gaussianity is consistently ac- 
counted for. We compute, for a generic barrier, the first two scale-independent halo 
bias parameters, the conditional mass function and the halo formation time probabil- 
ity, including the effects of non-Gaussianities. We also provide the expression for the 
two-constant-barrier first-crossing rate when non-Markovian effects are induced by a 
top-hat filter function in real space. 

Key words: cosmology: theory - large scale structure of the universe 



1 INTRODUCTION 

The distribution in mass of dark matter halos, as well as 
their clustering properties, formation history, and merging 
rate, play an important role in many problems of modern 
cosmology, because of their relevance to the formation and 
evolution of galaxies and clusters, and of their sensitivity to 
the statistical properties of the primordial density field. In 
particular, the most massive halos evolved from rare fiuc- 
tuations in the primordial density field, so their abundance 
and clustering properties are sensitive probes of primordial 
non-Gaussianities (Matarrese et al. 1986; Grinstein & Wise 
1986; Lucchin et al. 1988; Moscardini et al. 1991; Koyama 
et al. 1999; Matarrese et al. 2000; Robinson & Baker 2000; 
Robinson et al. 2000; LoVerde et al. 2008; Maggiore & Ri- 
otto 2010c; Lam & Sheth 2009; Giannantonio & Porciani 
2010), which could be detected or significantly constrained 
by various planned large-scale galaxy surveys, see, e.g. Dalai 
et al. (2008) and Carbone et al. (2008). Furthermore, the 



primordial non-Gaussianities (NG) alters the clustering of 
dark matter halos inducing a scale-dependent bias on large 
scales (Dalai et al. 2008; Matarrese & Verde 2008; Slosar 
et al. 2008; Afshordi & ToUey 2008) while even for small pri- 
mordial NG the evolution of perturbations on super-Hubble 
scales yields extra contributions on smaller scales (Bartolo 
et al. 2005; Matarrese & Verde 2009; Bartolo et al. 2010). 
The halo mass function can be written as 



dn{M) 
dM 



p dlncr~^(M) 
M2 dlnM 



(1) 



where n{M) is the number density of dark matter halos of 
mass M, a{M) is the variance of the linear density field 
smoothed on a scale R corresponding to a mass M, and p 
is the average density of the universe. The basic problem is 
therefore the computation of the function /(cr). Analytical 
computations of the halo mass function are typically based 
on Press-Schechter (PS) theory (Press & Schechter 1974) 
and its extension (Peacock & Heavens 1990; Bond et al. 
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1991) known as excursion set theory (see Zentner (2007) 
for a review). In excursion set theory the density perturba- 
tion depends stochastically with the smoothing scale, and 
the problem of computing the probability of halo formation 
is mapped into the so-called first-passage time problem in 
the presence of a barrier. With standard manipulations (see 
e.g. Zentner (2007)), the function /(cr) which appears in (1) 
is related to the first-crossing rate J" by f{a) — 2a^T{a^). 

In a recent series of papers (Maggiore & Riotto 
2010a,b,c) (hereafter MRl, MR2 and MRS, respectively), 
the original formulation of excursion set theory has been ex- 
tended to deal analytically with the non-Markovian effects 
which are induced either by the use of a realistic filter func- 
tion, or by non-Gaussianities in the primordial density field. 
In the original formulation of Bond et al. (1991) the prob- 
lem with the density field smoothed using a top-hat filter 
in wavenumber space was solved analytically and numeri- 
cal techniques wore adopted for the case of non-Markovian 
noises. The use of a top-hat window function in momen- 
tum space has the technical advantage that the evolution 
of the smoothed density field with the smoothing scale be- 
comes Markovian, but its important drawback is that is it 
is not possible to associate a well-defined mass to a region 
smoothed with such a filter (see Bond et al. (1991); Zentner 
(2007); Maggiore & Riotto (2010a)). For any other choice of 
filter function such as a top-hat function in real space (for 
which the relation between the mass M and the smooth- 
ing scale R is well-defined and is simply M = (4/3)7ri?''p) 
the actual evolution of the smoothed density field with R is 
non-Markovian. 

The same happens if the initial conditions for the grav- 
itational potential and/or the density contrast are non- 
Gaussian and the problem was solved first in MRS. The 
basic idea is to rcforumlate the first-passage time problem 
in the presence of a barrier in terms of the computation of a 
path integral with a boundary (i.e. over a sum over all "tra- 
jectories" S(S) that always stay below the barrier), and then 
to use standard results from quantum field theory and statis- 
tical mechanics to express this path integral in terms of the 
connected correlators of the theory. This allows us to include 
the effect of non-Markovianities arising, e.g., from the non- 
Gaussianities. In particular, in MRS we have shown how to 
include the effect of a non-vanishing bispectrum, while the 
case of a non-vanishing trispectrum was considered in Mag- 
giore & Riotto (2010d) (see also D'Amico et al. (2010) for 
an approach to non-Gaussianities which combines our tech- 
nique with the saddle point method developed in Matarrese 
et al. (2000)). 

An essential ingredient of excursion set theory is a 
model for the collapse of a dark matter halo. In its simplest 
implementation, one uses the spherical collapse model. This 
model, however, is certainly a significant over-simplification 
of the complicated dynamics leading to halo formation and 
can be improved in different, complementary, ways. A cru- 
cial step was taken by Sheth, Mo & Tormen (2001) who took 
into account the fact that actual halos arc triaxial (Bardeen 
et al. 1986; Bond & Myers 1996) and showed that an ellip- 
soidal collapse model can be implemented, within the excur- 
sion set theory framework, by computing the first-crossing 
rate in the presence of a barrier Bst(S), 



1 + OA 



(2) 



which depends on 5* = cr'^ ("moving barrier"), rather than 
taking the value 5c{z) of the spherical collapse, which is red- 
shift dependent, but independent of S. Physically this re- 
fiocts the fact that low-mass halos (which corresponds to 
large S) have larger deviations from sphericity and signifi- 
cant shear, that opposes collapse. 

Notice that, to improve the agreement between the pre- 
diction from the excursion set theory with an ellipsoidal col- 
lapse and the N-body siirmlations, Sheth, Mo & Tormen 
(2001) also found that it was necessary to multiply Sc{z) by 
^/a, where ^/a ~ 0.84 was obtained by requiring that their 
mass function fits the GIF simulation. In MR2 we proposed 
a physical justification for the introduction of this parame- 
ter in the halo mass function, suggesting that some of the 
physical complications inherent to a realistic description of 
halo formation could be included in the excursion set the- 
ory framework, at least at an effective level, by treating the 
critical threshold for collapse as a stochastic variable, whose 
scatter reflects a number of complicated aspects of the un- 
derlying dynamics (see also Audit et. al. (1997); Lee & Shan- 
darin (1998); Sheth, Mo & Tormen (2001) for earlier related 
ideas). Solving the first-passage time problem in the pres- 
ence of a barrier which is diffusing around its mean value, 
it was found in MR2 that the coefficient a can be related 
to the diffusion coefficient Db oi the stochastic barrier as 
a = 1/(1 + Db)- The numerical value of Db, and therefore 
the corresponding value of a, depends among other things 
on the algorithm used for identifying halos. From recent N- 
body simulations that studied the properties of the collapse 
barrier, a value Db — 0.25 was deduced in MR2, predicting 
a ~ 0.80, in agreement within the accuracy of the MRl pre- 
diction (~ 20%) with the value of a extracted directly from 
a fit to the mass function (see also Corasaniti & Achitouv 
(2010) for recent related work). 

The path-integral formulation developed in MRl and 
MRS WEis restricted to the case of a constant barrier 5c{z) 
and it was subsequently generalized to the case of the el- 
lipsoidal moving barrier in De Simone et al. (2010). In the 
present paper we further develop the path-integral formula- 
tion of the excursion set theory to calculate, for a generic 
moving barrier and for Gaussian and non-Gaussian initial 
conditions, other basic quantities necessary to characterize 
the physics of dark matter halos like halo bias, accretion 
rates, formation times, merging, halo assembly bias and so 
on. 

We know that dark matter halos typically form at sites 
of high density peaks. The spatial distribution of dark mat- 
ter halos is therefore a biased tracer of the underlying mass 
distribution. A standard way to quantify this difference be- 
tween halos and mass is to use a bias parameter bh, which 
can be defined as the ratio of the ovordensity of halos to 
mass, or as the square root of the ratio of the two-point 
correlation function (or power spectrum) of halos to mass. 
Like the halo mass function, analytic expressions for the halo 
bias can bo obtained from the excursion set theory based on 
the spherical gravitational collapse model (Cole & Kaiser 
1989; Bond et al. 1991; Mo & White 1996) and for the el- 
lipsiodal one (Sheth, Mo & Tormen 2001). The approach 
to the clustering evolution is based on a generalization of 
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the so-called peak-background split (Bardeen at al. (1986)) 
which basically consists in splitting the mass perturbations 
in a fine-grained (peak) component filtered on a scale R 
and a coarse-grained (background) component filtered on a 
scale Ro S> R- The underlying idea is to ascribe the col- 
lapse of objects on small scales to the high frequency modes 
of the density fields, while the action of large-scale struc- 
tures of these non-linear condensations is due to a shift of 
the local background density. In the excursion set theory 
the problem of computing the probability of halo formation 
is mapped into the first-passage time problem of a random 
walk which starts from a given value of the density contrast 
So at a given radius Rq corresponding to a given value of the 
variance a{Mo). When the random walk performed by the 
smoothed density contrast is Markovian, the first-crossing 
rate is easily computed by a simple shift of the initial con- 
ditions. This is due to the fact that, being the noise white, 
the memory about the way the system arrived at the point 
So at a given time is lost. On the contrary, when the random 
walk is non-Markovian, the system has memory effects and 
it remembers how it arrived at So- This influences the subse- 
quent first-crossing rate. The computation of the halo bias 
mass function in the case in which the non-Markovianity 
is induced by the choice of a top-hat window function in 
real space, and within a spherical collapse model, has been 
recently performed in Ma et al. (2010). In this paper we 
perform the calculation of the halo bias parameters for the 
ellipsoidal barrier and when non-Gaussian initial conditions 
introduce non-Markovianity, see also Giannantonio & Por- 
ciani (2010) for a treatment of NG halo bias. 

The excursion set theory is also a powerful formalism 
for studying the formation history of halos. The most im- 
mediate quantity of interest is the conditional mass func- 
tion. Given a halo of mass Mo at redshift Za, one can com- 
pute the average manner in which this mass was partitioned 
among smaller halos at some higher redshift Zb > Za- The 
conditional mass function is simply the average number 
of halos of mass M„ at redshift zi, that are incorporated 
into an object of mass Mo at redshift Za- In the language 
of excursion set theory this can be formulated as a two- 
barrier problem, i.e. in terms of the conditional first cross- 
ing rate, T{Bb{Sn), Sn\Ba{So), So), describing the rate at 
which trajectories make their first crossing of the barrier 
Bi,{S) = B{S,z = Zb) at a value S = Sn, corresponding to 
the mass M„, under the condition that, at an earlier "time" 
S — So (corresponding to the mass Mo; recall that decreas- 
ing the variance S the corresponding mass Af (5) increases, 
so So < Sn means Mo > M„), they crossed the threshold 
Ba{S) = B{S,z = Za). Then, a halo of mass Mo has its 
mass partitioned on average among a spectrum of halos at 
redshift Zb as (Lacey & Cole 1993; Zentner 2007) 



dn{Mn\Mo) _ Mo 



dMn 



-Tf- -^(BbiSn), Sn\Ba{So), So) 



dSn 



dM„ 



(3) 



The function F{Sn, Bb{Sn)\So, Ba{So)) gives the probability 
of the second barrier first-crossing at a particular value of 
Sn, while the factor (Mo/Mn) converts it from a probability 
per unit mass of halo Mo into the number of halos of mass 
Mn- The two- barrier result can also be manipulated to yield 
the average mass accretion rate, halo formation time and so 
on. 

The relationship between the unconditional mass func- 



tion and the first-crossing distribution associated with 
barrier-crossing random walks has been extended to obtain 
the conditional mass function of halos by Bond et al. (1991) 
and Lacey & Cole (1993) within the spherical collapse (the 
so-called extended Press-Schechter model). The two- barrier 
first-rate probability has a simple analytic form in the con- 
stant barrier spherical collapse model. Again, this is because 
the random walk performed by the smoothed density con- 
trast is Markovian. For a moving barrier (such as the ellip- 
soidal collapse model), however, exact analytic forms have 
been found only for the special case of a linear barrier (Sheth 
1998) while Lam & Sheth (2009) have proposed a better mo- 
tivated approach to the problem by offering a Taylor series- 
like approximation for a general moving barrier, see also 
Giocoli et al. (2007). Zhang et al. (2008) have provided 
analytical expressions for the two-barrier first-crossing rate 
for the ellipsoidal collapse and Gaussian initial conditions. 
In this paper we compute this conditional probability using 
the path-integral formulation for a generic moving barrier for 
Gaussian and non-Gaussian initial conditions. A by-product 
of such a calculation is the determination of the halo forma- 
tion time probability. 

The paper is organized as follows. In section 2 we sum- 
marize the basic ingredients for the calculation of the first- 
crossing rate from the excursion set theory and a generic 
moving barrier. In section 3 we compute the conditional 
probability necessary to deduce the halo bias parameters, 
both in the Gaussian and non-Gaussian case. Section 4 con- 
tains the computation of the two-barrier first-crossing rate 
for a generic moving barrier and again for both Gaussian 
and non-Gaussian initial conditions. In section 5 we present 
our results for the halo formation time probability. Finally, 
Section 6 contains our conclusions and a summary of the 
main results, while some technical material is collected in 
the Appendices. In particular. Appendix A contains some 
useful numerical fits, while Appendix B contains the com- 
putation of the two-barrier first crossing rate including the 
non-Markovian effects coming from the choice of a top-hat 
filter in real space. 



2 PATH INTEGRAL FORMULATION OF 
EXCURSION SET THEORY FOR A 
MOVING BARRIER 

Let us discuss the basic points of the original fornmlation 
of excursion set theory for a moving barrier. We will closely 
follow MRl and De Simone et al. (2010); at the expense of 
being ripetitive, we will report here various details that the 
reader can find in these references. This will hopefully help 
to follow and speed up the calculations of the subsequent 
sections. 

In the excursion set theory, one considers the density 
field S smoothed over a radius R, and studies its stochastic 
evolution as a function of the smoothing scale R. As it was 
found in the classical paper by Bond et al. (1991), when the 
density S{R) is smoothed with a sharp filter in wavenumber 
space, and the density fluctuations have Gaussian statistics, 
the smoothed density fleld satisfles the equation 



dSjS) 
dS 



■ viS) , 



(4) 
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where S = <y^{R) is the variance of the hnear density field 
smoothed on the scale R and computed with a sharp filter in 
wavenumber space, while ri{S) is a stochastic variable that 
satisfies 



(5) 



where 5d denotes the Dirac delta function. Equations (4) 
and (5) are the same as a Langevin equation with a Dirac- 

dclta noise r](S), with the variance S formally playing the 
role of time. Let us denote by n((5, S)d5 the probability den- 
sity that the variable 5{S) reaches a value between 5 and 
5 + d5 hy "time" S. In the general non-Markovian case it 
is not possible to derive a simple, local, differential equa- 
tion for n((5, S) (indeed, it can be shown that \l(5, S) rather 
satisfies a complicated integro-differential equation which is 
non-local with respect to "time" S, see eq. (83) of MRl), 
so one cannot proceed as in the Markovian case where, as 
we will review below, n((5, S) is determined by the solution 
of the Fokker-Planck equation with appropriate boundary 
conditions. Rather, we construct the probability distribu- 
tion n((J, S) directly by summing over all paths that never 
exceeded the corresponding threshold, i.e. by writing Yl{5, S) 
as a path integral with a boundary. To obtain such a repre- 
sentation, we consider an ensemble of trajectories all starting 
at So = from an initial position 5{0) = Sq and we follow 
them for a "time" S. We discretize the interval [0, S] in steps 
AS = e, so Sk ~ ke with k = 1, . . .n, and Sn = S . A trajec- 
tory is then defined by the collection of values {(5i, . . . , (5„}, 
such that 5{Sk) = Sk and B{Si) = Bi. The probability den- 
sity in the space of trajectories is 

W{5,u Si, J„; 5n) = {Sd{S{Si)-Sj) . . . <5d(<5(S„)-<5„)) ,(6) 

where Sd denotes the Dirac delta. Then the probability of 
arriving in Sn in a "time" Sn, starting from an initial value 
So, without ever going above the threshold, is (Bond et al. 

(1991)) 

rBi i-B„ 
n^biS„;S„) = I dSi... dSn-1 (7) 



dSi . . . / dSn- 

-oo J —oo 

X W^((5o; . . . , 5„-i, 5„; Sn). 



The label "mb" in limb stands for moving barrier. The func- 
tion W{So; (5i, . . . , Sn-i,Sn; Sn) can be expressed in terms of 
the connected correlators of the theory. 



W{5o;5x,...,5n;Sn) = j VXe^ , 
where 

/'_,,_ f°° d\i dXn 

J J — OO 

and 

Z = i ^ \iSi 



(8) 



(9) 



(10) 



n n 



p=2 



P 



il = l ip=l 



Here {5i...Sn}c denotes the connected n-point correlator. 
So 



When S{S) satisfies eqs. (4) and (5) (which is the case 
for sharp filter in wavenumber space) the two-point function 
can be easily computed, and is given by 



(S{Si)S{Sj))=mm{Si,Sj). 



(12) 



In the rest of this section we will restrict ourselves to the 
Gaussian and Markovian case. 



2.1 The case of Markovian noise 

Taking the derivative with respect to the time S„ = S 
of eq. (11) and using the fact that, when multiplying 
exp{i Ai(5i}, iXj {j = l,---,n) can be replaced dj = 
d/dS^, we discover that IImb(<5n;5„) satisfies the Fokker- 
Planck (FP) equation 



anmb(5n;S„) ^ 19'nmb(<5„;5'„) 
2" 



dSn 



dS^ 



(13) 



In the continuum limit, the boundary condition to be im- 
posed on the solution of eq. (13) is 



ninb(<5„; Sn) = for S„ ^ B„ 



(14) 



In the continuum limit the first-crossing rate is then given 
by 



~ dS 

J —a 



ddn limhi^n; Sn) 



f 

o —a 



(15) 

dSn 



The first term on the right-hand side vanishes because of the 
boundary condition, while the second term can be written 
in a more convenient form using the FP equation (13), so 



mh{Sn) — 



dS 



d^n^y,{Sn;Sn) 



dS^ 



1 Snmb(<5„; Sn) 
"2 



dSn 



(16) 



To compute the probability Ylmh{Sn, Sn) we proceed in the 
following way. At every i-th stop of the path integral we 
Taylor expand the barrier around its final value, £is first 
suggested by Lam & Sheth (2009), 



Bi 



Sn)^ , 



(17) 



where 
= 



(p) _ d^B{Sn) 



dSl ' 

(so in particular Bn'' = B{Sn))- We now perform a shift in 
the integration variables Si {i 
integral 



(18) 
ft in 

1, . . . , n — 1) in the path 



,(p) 



Si^Si-J2^ (Si- Sn)" 



Then UmbiSn', Sn) can bo written as 



/Bl ("Sn-l r pBn pBn /• 

dSi... j dSn-i / OAe^. (11) n^h{Sn;Sn)= dSl... dSn-l / 



(19) 



(20) 
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where 



n 1 ^ 

Z = i^^Ai(5i — - mm{Si, Sj 



i=l i,j=\ 

Jz3 ^ r(p) 



(21) 



We next expand 

{ n — 1 oc j^i.P) 

exp < i ^ A^ ^ — ^ - 5*^)^ 



n — 1 oo 



i— 1 p— 1 



(22) 



-- ^i^J 2^ ^j^i {Si - Sn)^ {Sj - Sn)" + ' 

and we write Ilinh{Sn', Sn) as 



(23) 



n^b{Sn;S„) = n^^(5„;5„)+n«(,5„;5„) 

+n^2{Sn;Sn) + --- . 

For the zero-th order term we can immediately take the 
continuum limit, using the results of MRl, and we get the 
standard probabihty density of excursion set theory in the 

Markovian and Gaussian case, 



nL°^(5„;S„) 



\/2'KSn 



(24) 



The terms and are given by 



nL'^(<5n;5„) 



and 



nL'^(5n;5„) 

X {Si - S^Y {Sj - S^Y didjW^-^{5o; di,...,6„; Sn) , (26) 



E/ 'd5,...dSn-.J2^ (25) 

i=l p=l 

X {Si - SnY diW^'^{5o; 5i, . . . , S„) , 



^ n-l .s„ oo 

= - ^ / d<5i...d<5„-i ^ 

^ «/ — OO 1 



where 

W^'^{5o;Si,...,dn;Sn) = 



(27re)"/2 



(27) 



and superscript "gm" (Gaussian-Markovian) reminds us 
that this value of W is computed for Gaussian fluctuations, 
and when the evolution with respect to the smoothing scale 
is Markovian. Their continuum limit is more subtle, and can 
be computed using the technique developed in MRl, as we 
review below. 

We have therefore formally expanded Ilinb{Sn, Sn) in a 
series of terms 11^^, 11^^, etc., in which each term is itself 
given by an infinite sum over indices p,q, . . . . We have to 
evaluate the continuum limit of objects such as 

rB„ 

F{Si) / dSi... d5n-i diW^{5o; 5i, . . . , 5„; Sn) , (28) 

■ 1 J —oo 



where F denotes a generic function. To compute this expres- 
sion we integrate di by parts, 

"Sr. 

dSi . ..d6n-i diW^{6o;Si, . . .,6n;Sn) 



I 

J — ( 



= / d5i...d5i...d5„-i (29) 

J —oo 

xW^((5o;<5i, ■ ■ ■ ,5i = Bn, ■ ■ .,5n-l, Sn, Sn) , 

where the notation dSi means that we must omit dSi from 
the list of integration variables. We next observe that VF^™ 
satisfies 

W^'^{6o;di, . . . ,6i = Bn, . . . ,Sn;Sn) 

= W^'^{S„:S-,,...,5,-j,Bn;Si) 

X W°'"(S„; S,+i, ...,5„;S„~ S^) , (30) 

as can be verified directly from its explicit expression (27). 
Then 

d5i...dSi^i / dSi+i . . .d6„-i 



(31) 



xW<''^{5o;di,...,5i-i,Bn;Si) 

XW^'^{Bn;di+l,...,5n;Sn-Si) 

= n«'"(5o; Bn, Si)n«'"(B„; Sn, Sn - Si) , 



and to compute the expression given in eq. (28) we must 
compute objects such as 



J2 F{Si)n^'^{5o; Bn, Si)U^{Bn; 5„; 5„ - Si) 



(32) 



We then need to know n^™((5o; B„; Si). By definition, in the 
continuum limit this quantity vanishes, since its second ar- 
gument is equal to the the threshold value S„. However, in 
the continuum limit the sum over i becomes 1/e times an 
integral over an intermediate time variable Si, 



i=l ^ 



(33) 



so we need to know how 11^'^ {So \ Bn; Si) approaches zero 
when e — ^ 0. In MRl it was proven that it vanishes as ^/e, 
and that 



n^-{So;Bn;Sn) = V~e ^^LZ^e-^"--'"^"/^''"^ 



\Yk S] 



i/2 



+ 0{e) 



(34) 



Similarly, for Sn < Bn, 

'/t^ Sn 

(35) 

In the following, we will also need the expression for 11^™ 
with the first and second argument both equal to Bn , which 
is given by (see again MRl) 



nS-(B„;B„;5) = 



2^S'J' 



(36) 



In order to finalize the computation, we must either perform 
some approximation, or identify a suitable small parameter, 
and organize the terms in a systematic expansion in such 
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a small parameter. In Dc Simone et al. (2010) we have 
discussed in detail two different expansion techniques (one 
based on a systematic expansion in derivatives for a slowly 
varying barrier, and the other in which a large number of 
terms are resummed), which were shown to provide very 
close numerical results. Furthermore, it was found that the 
results obtained with these systematic expansions are in the 
end numerically very close to that obtained with a simpler 
albeit more empirical procedure, which amounts to approxi- 
mating {S„ — SiY'^ ~ S^~^ inside the integrals in eqs. (25) 
and (26), and at the same time truncating the sum over p 
in eq. (25) to p = 5 (while, in this approximation, n|^^ does 
not contribute). This is in fact equivalent to the approxima- 
tion made by Sheth & Tormcn (2002) (ST in the following) 
and gives the same results. However, as discussed in Sec- 
tion 3.1 of De Simone et al. (2010), when one makes the 
approximation [Sn — Si)''^^ ~ S^^^, one must also neces- 
sarily truncate the sum to a maximum value, otherwise the 
first-crossing rate resums to a trivial result, where all cor- 
rections due to the ellipsoidal barrier disappear. Therefore, 
the procedure of replacing {S„ — Si)''~^ — > SZ~^ inside the 
integrals and, at the same time, truncating the sum, must 
be viewed as a simple heuristic procedure to get a result 
which is numerically close to the result of more systematic 
expansions. Since this procedure is technically much simpler 
than the systematic expansions discussed in De Simone et 
al. (2010), and works well numerically, we will adopt it in 
the following. 

We first compute H^^. Before performing the above ap- 
proximation, the expression of n^^(5„;S'„) in eq. (25) can 
be rewritten as 



^mh) terms which arc neglected give contributions pro- 
portional to higher powers of (i?„ — 5n)- Since in the end the 
mass function is obtaned from the first-crossing rate (16), we 
actually only need the first derivative of TlmhiSn; Sn) eval- 
uated at d„ = Bn, and terms proportional to {B„ — 5„)'^ 
with iV > 2 give a vanishing contribution. 

Higher-order contributions to the first-crossing rate 
vanish. In fact, in the same approximation one finds (De 
Simone et al. 2010) that H^,^^^' vanishes as (S„ — (5„)" for 
S„ Bn, so its first derivative dH^^^'^^ /dSn evaluated in 
Sn = Bn, which according to eq. (16) gives its contribution 
to the first-crossing rate, vanishes for all n ^ 2. 

The total first-crossing rate for a moving barrier, in the 
approximation discussed above, is therefore given by 



^ mb(5'n) — 



2nSl 



.3/2 



E 

p=0 



-S^/(2S„) 
1^ 



■{Bn.+V{Su)) , 



(40) 



where 



V{Sn) 



E 



{-SnY d^Bn 

v\ dSl 



(41) 



When applied to the ellipsoidal barrier given in eq. (2), one 
recovers the ellipsoidal collapse result of Sheth & Tormen 
(2002). 



nL'^(5n;S„) = 



Bn{Bn-6u) (-IF 



— On) 



B. 



(p) 



(Sn - s,r-^'/'^ 



Jo 



xe 



-Bl/(2Si) -{Bn-6nr/[2{Sn-Si)] 



(37) 



Since this integral is finite in the limit Sn — > Bn, taking the 
approximation {Sn — SiY~^ ~ {SnY^^ does not alter the 
convergence properties of the integral, but simplifies signif- 
icantly its computation, since 



dSi 



„ -Bi/{2Si) -(B. 

xe "'^ "e ^ 



V2^ 1 



Sf^{Sn - SiY/^ 

S„f/C2(S„-Si)) 

{2Bn - dnf ' 



- Bn ^y^-^^P 1^ 2Sn 

SO in this approximation H^^((5„; Sn) is given by 

2{Bn - Sn) ^-(2Bn-S„f/(2Sn) 



f {2Bn - ^n)' ] 
\ 2Sn J 



(38) 



2nS'i^^ 



E 

p=i 



>(p) 



(39) 



where the superscript "ST" reminds us that we have per- 
formed the approximations that are equivalent to those 
which give the ST mass function. A reason why this ap- 
proximation works well is that (at least for what concerns 



J^st{S„) 



aS^{z) ^-Bl/(2Sn) 



2i,sT 



l-h 



+0.4 ^(-1)^ 

faSc{z) ^-bI/(2S) 



0.6 



Sn 



aSc{z) 



2-kS\ 



.3/2 



1 + 0.067 



Sn 

a5l{z) 



(42) 



As it is well-known (Sheth & Tormen (2002)), this first- 
crossing rate is not normalized to unity. This is a basic dif- 
ference between the moving barrier and the constant (spher- 
ical) barrier model. When the barrier height is constant, all 
random walks arc guaranteed to cross the barrier because 
the rms height of random walks at Sn is proportional to 
\fS^. At sufficiently large Sn , all walks will have crossed the 
constant barrier. In the moving barrier case, in which the 
barrier diverges when Sn — ^ oo, not all trajectories intersect 
it. This is because the rms height of the random walk grows 
more slowly than the rate at which the barrier height in- 
creases and there is no guarantee that all random walks will 
intercept the barrier. It seems reasonable to associate the 
fraction of random walks that do not cross the barrier with 
the particles that in N-body simulations are not associated 
to bound states (Sheth & Tormen (2002)). 

After this rather long and technical summary of how to 
compute the first-crossing rate for a generic moving barrier, 
we are ready to compute conditional probabilities. 
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3 HALO BIAS 

We now apply the technique of the previous section to the 
computation of the halo bias, including the non-Markovian 
corrections coming from the NG. We will use a top-hat win- 
dow function in wavenumbcr space. The calculation of the 
non-Markovian effects on the bias from a top-hat window 
function in real space can be found in Ma et al. (2010). 



3.1 Conditional probability: the moving barrier 
case and Gaussian initial conditions 

We begin our analysis with the simpler case in which the 
density field is Gaussian. Since we are also taking a top- 
hat filter in wavenumber space, the evolution is the Marko- 
vian. To compute the bias, we need the probability of form- 
ing a halo of mass M, corresponding to a smoothing ra- 
dius R, under the condition that the smoothed density 
contrast on a much larger scale Rm has a specified value 
Sm = S{Rm)- We use J^inh{Sn\Sm, Sm) to denote the cor- 
responding conditional first-crossing rate, i.e., the rate at 
which trajectories first cross the barrier 5 = B{S) at time 
S„, under the condition that they passed through the point 
5 = Sm at an earlier time Sm- We also use the notation 

J^inh{S„\0) = J^inh{S„\5m = 0, Sm = 0), SO I'mbiSnlO) is the 

first-crossing rate when the density approaches the cosmic 
mean value on very large scales. 

The halo overdensity in Lagrangian space is given by 
(Mo & White (1996); see also Zentner 2007 for a review) 



J- + Ohalo ■ 



:^mb(s„|o) 



(43) 



The relevant quantity for our purposes is the halo condi- 
tional probability 

nhalo(^n, 'S'nl^m, Sm^ = 

dSi...d6m... f""-' dSn-iW (<5o = 0; <5i, . . . , 5„; S^) 
Clo dSm-iW {So = 0; 5i, • • • , Sm] Sm) 

(44) 

where the hat over dSm means that dSm must be omitted 
from the list of integration variables. The numerator is a 
sum over all trajectories that start from So — at S = 0, 
have a given fixed value Sm at Sm, and a value Sn at Sn, while 
all other points of the trajectory. Si, ... , Sm-i,Sm+i, • • • Sn-i 
are integrated up to the corresponding value of barrier, and 
we use the notation Bi = B{Si). The denominator gives the 
appropriate normalization to the conditional probability. 

The conditional first-crossing rate J^inb{Sn\Sm, Sm) 
is obtained from the conditional probability 

^hB.\o{Sn,Sn\Sm,Sm) USiug 

d f"" 

J^inh{Sn\Sm, Sm) = —-^^ / dSn^haloi^n, Sn\Sm, Sm) -(^S) 
OOn J_^ 

Since we axe considering the Gaussian case, with a top-hat 
filter in wavenumber space, the probability density W fac- 
torizes. 



W^'^{5o-Sl,...,5m,...,5n;Sn) 
W'''^{So;Sl,...,5m-l,Sm;Sm) 
^W^'^iSm-, 5m + l, ...,Sn;Sn- Sm) , 



and the halo probability nhaio(5n, Sn\Sm- Sm) in eq. (44) be- 
comes identical to the probability of arriving in Sn at time 
Sn, starting from Sm at time Sm for the moving barrier, 
reflecting the fact that the evolution of S{S) is in this case 
Markovian. We can then compute Ilhaio as in the previous 
section, performing a shift of the remaining integration vari- 
ables Si with i = (m -h 1, • • • , n — 1), 



Si — >■ S, 



_ \ 



{Si — SnY^ 1 



(47) 



and we get 



-(5„-«™)V{2{S„-S„)) 



•\/27r(S'„ — Sm) ^ 

-(2S„-5„-5,„)V(2(S„-S,„))> 



+ 



2{Bn - Sn) 



\r2^{Sn - Sm)'"-'^ 

2{Bn - Sn)' 
^{Sn-Sm)^'^' 



-(2Br,-6n-Smr/C2(S„-Sm))T> 
" I mn 



-(2S„ 



-«„)V(2(S„-S„)).p2 

' TTi 



where 

Vmn = ViSm, Sn) = V (^HLZ^B^) 



(48) 



(49) 



p=i 



In the following, we will use this quantity with the sum trun- 
cated to p = 5, as discussed in Section 2 and as advocated by 

Sheth & Tormen (2002) for the conditional mass function. 

The calculation of the conditional first-crossing rate 
J^inh{Sn\Sm, Sm) proceeds by taking the derivative with re- 
spect to Sn 



(Bn 



+ ^(^-^")e-^^.(50) 



^mbiSr.\Sm,Sm)- ^^^^ _ ^^^^^^ 

In a sufficiently large region Sm <S Sn and Sm Sn. Then, 

expanding to quadratic order in Sm and after mapping to 
Eulerian space, we find the first two Eulerian bias coefficients 

Bn 1 



and 



62 



B. 



Sn B„+V{0,Sn) 



Bn 



Sn Sn iB„+V{0,Sn))Sn 



(51) 



(52) 



The above results hold for a generic barrier. We now examine 
it for different collapse models. 

3.1.1 Constant barrier 

We first apply these results to the spherical collapse model 

where the barrier is constant Sc{z). Using the standard no- 
tation f = Sc{z)/a we get, for the bias coefficients, 

2 1 

(53) 



(46) 



i-'C 

and 



(54) 
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Figure 1. The halo bias at linear order as a function oiv = 5c/cr- 
Our formula in Eq. (55) (solid black line) is compared to other 
results: the fit to N-body simulations as in Tinker et al. (2010) 
(dashed red), and two standard approximations: spherical (dotted 
blue) and the SMT expression for the bias (dot-dashed magenta), 
which is also predicted by the diffusing barrier model of Section 
3.1.3. 



3.1.2 Ellipsoidal harrier 

For the the ellipsoidal barrier (2), we get, for the bias coef- 
ficients. 



6i 

and 

h2 



H-0.4 
1 



f-V 



1 + 0.067 



' 1 



(55) 



'51 

..2 



1 + 0.4 



f-V 



V 

5f 



1 0.067 { ^ 



(56) 



For large masses, 



1, the bias coefficient scales like 



(57) 



6i ~ ^ — . 

Oc 

Observe that this result differs from that proposed by Sheth, 
Mo & Tormen (2001), which in the large mass limit rather 
scales like h\ ~ a{v^ /Sc), i.e. it is smaller than (57) by a 
factor y/a. The proportionality to ^/a in eq. (57) can be 
traced to the fact that, in the large mass limit, bi is dom- 
inated by the term B„/Sn in eq. (51) and, for the barrier 
given in eq. (2), B„ is proportional to y/a. In other words, 
this dependence is a consequence of the fact that, in the 
large mass limit, the barrier of Sheth, Mo & Tormen (2001) 
does not reduce to the spherical collapse model barrier Sc, 
but rather to ^/aSc, so the large mass limit (57) can be for- 
mally obtained from the spherical collapse model by per- 
forming the replacement 5c — > \/a5c. For Sheth, Mo & Tor- 
men (2001), the physics of spherical collapse always gives 
(5c, not y/aSc, so it is really only the variances which were 
rescaled, — >■ cP' ja. Since the unconditional mass func- 
tions are the same in our case and in the one discussed by 
Sheth, Mo & Tormen (2001), but the bias factors are not, 
this means that one could, in principle, discriminate between 
whether a should be thought of as rescaling or cP' . 



In Fig. 1 we show our result (55) for the bias and com- 
pare it to the fit to N-body simulations of Tinker et al. 
(2010) and to two standard approximations: spherical col- 
lapse and the SMT result. Our result for the bias is mani- 
festly in better agreement with N-body simulations than the 
commonly-used approximations. 

3.1.3 The diffusing barrier 

In MR2 it has been proposed a model in which the barrier 
performs a diffusing motion, with diffusion coefficient Db, 
around an average value which for the spherical collapse 
model is simply 5c{z) (without any factor of y'Ti) while for 
the ellipsoidal model is given 



1 + 0.4 



(58) 



i.e. by eq. (2) without the factor ^/a. Equation (58) is in fact 
the barrier which actually follows for an ellipsoidal model for 
collapse, while the factor in eq. (2) was simply inserted 
by hand, in order to fit the data. Including the stochastic 
motion of the barrier, which is meant to mimic a number of 
random effects in the process of halo formation, it was found 
in MR2 that the relative motion of the trajectory S{S) and 
of the barrier is a stochastic motion with diffusion constant 
{1 + Db) so that, for instance, in the Markovian case the 
evolution of the probability distribution is governed by a 
Fokker-Planck equation of the form 



9n^b(5n; S„) _ (1 + Db) 9'n^b(<5n; S„ 



dS„ 



(59) 



The result for this diffusive barrier can therefore be obtained 
from the result for the spherical collapse barrier, or from the 
result for the barrier (58), by formally rescaling 5* — >■ {l/a)S, 
where a = 1/(1 + Db). For the halo mass function, which 
depends only on the combination ly — Sc/a — 5c /^fS, the 
rescaling Sn — >■ {l/a)Sn is equivalent to the rescaling 5c — > 
^/a5c, and therefore the diffusive barrier model produces the 
same result as that obtained with the SMT barrier (2). 

This equivalence, however, does not extend to the halo 
bias. In fact, in the large mass limit, the halo bias deduced 
using the barrier (58) in eq. (51) is foi ~ {^^/5c) = 5c/Sr., 
and rescaling Sn — >■ {l/a)Sn gives 



fei ~ a — 



(60) 



which differs from eq. (57) by a factor ^/a and agrees with 
Sheth, Mo & Tormen (2001). 



3.2 Adding non-Markovianity from a top-hat 
window function in real space, in the high 
mass limit 

Furthermore, one should add to this result the non- 
Markovian corrections due to the use of a top-hat filter 
function in coordinate space. In the large mass limit, where 
eq. (58) reduces the constant barrier, the computation of 
the bias was performed in Ma et al. (2010), where it was 
found that bi gets multiplied by a factor 1/(1 — a«:), so for 
the diffusive barrier model with Markovian corrections due 
to the filter we get, in the large mass limit. 
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61 



1 V 
an 6c 



1 — an 



^/a ■ 



(61) 



where /t is a parameter that controls the non-Markovian 
effects of the filter function, whose numerical value can be 
estimated as in MRl, but which for accurate fitting is better 
treated as a free parameter, see Ma et al. (2010). Observe 
that in Ma et al. (2010) the spherical collapse result was 
rescaled according to 5c — > '\/aSc, which actually corresponds 
to the replacement in the SMT barrier model and not to the 
diffusive barrier model, so the resulting value for the bias 
was incorrectly obtained to be bi ~ [^/{l — aft)] (i^^ /(5c)- 

In particular, consider that the rescaling k — > a/t, which 
leads to the factor (1 — an) in the denominator of eq. (61), is 
obtained assuming that the barrier performs a simple diffu- 
sive Markovian motion. In a realistic description, the actual 
stochastic motion can be more complicated and the reader 
should be aware of the limitations of this procedure. Unfor- 
tunately, we have not been able to find an analogue of the 
expression (61) for small masses. 



3.3 Conditional probability: the moving barrier 
case and non- Gaussian initial conditions 

Deviations from Gaussianity are encoded, e.g., in the con- 
nected three- and four-point correlation functions which are 
dubbed the bispectrum and the trispectrum, respectively. 
A phenomenological way of parametrizing the level of NG 
is to expand the fully non-linear primordial Bardeen grav- 
itational potential $ in powers of the linear gravitational 
potential $l 



$ = $L + /nL - 



(62) 



The dimensionless quantity /nl sets the magnitude of the 
three-point correlation function (Bartolo et al. (2004)). If 
the process generating the primordial NG is local in space, 
the parameter /nl in Fourier space is independent of the 
momenta entering the corresponding correlation functions; 
if instead the process which generates the primordial cos- 
mological perturbations is non-local in space, like in models 
of inflation with non-canonical kinetic terms, /nl acquires 
a dependence on the momenta. The strongest current limits 
on the strength of local NG set the /nl parameter to be in 
the range —4 < /nl < 80 at 95% confidence level (Smith et 
al 2010). The goal of this subsection is to compute the halo 
bias parameters in the presence of NG and for the ellipsoidal 
collapse, using the technique developed in MRS, which gen- 
eralizes excursion sot theory to deal with non-Gaussianitics 
in the primordial density field. We will also work with a 
top-hat window function in real space. 

Similarly to the Gaussian case, the probability of ar- 
riving in 6n in a "time" S„, starting from the initial value 
6o = 0, without ever going above the threshold, in the pres- 
ence of NG is given by 

dSi . . . d6n-i 

-oo J —oo 

xWng(5o;5i, . . . ,5n-i,5„;S„). (63) 



where 



Wng('5( 



o; <5i, . . . , 5n; Sn) = j T>\ 



In 1 ^ 

: exp < i XiSi — - AjAj min(S'f , Sj 



X exp ■ 



(64) 



i,i,k=l 



and we have retained only the three-point connected cor- 
relator (5idjdk}c as a signal of NG. It is now clear where 
the non-Markovianity is coming from when non-Gaussian 
initial conditions are present: expanding H^ng in powers of 
j ,.^i{SiSj5k)c^i^j^k and going to the continuum, one 
obtains integrals over the intermediate times which intro- 
duce memory effects. 

As in eq. (44), the conditional probability relevant to 
the halo bias is obtained keeping Sm fixed, rather than 
treating it as an integration variables. In principle the non- 
Maxkovian contribution to the halo probability from NG 
should be computed separating the various contributions to 
the sum according to whether an index is smaller than or 
equal to m, larger than m and smaller than or equal to n. 
Fortunately, to compute the halo bias parameters, at the 
end one needs to take the limit Sm >C Sn and 5m <S 5n- 
In this limit we can the safely neglect the contribution to 
the sum from all indices running from 1 to m, since in 
general (5^5^) scales like (Sm^S?,,^^), where p,q ^ and 
p+q = 3, and therefore vanishes for Sm — >■ 0. In other words, 

in eq. (64) we can replace E"i,fc=i ^^^^ EL,fc=m+r ^s 
usual, we then expand the exponential and use the fact that 
i\k exp{i XiSi} = dk exp{i \i5i}^ and we use the fac- 
torization property (30) (Maggiore & Riotto (2010c)). This 
gives 



Wng((5o; 5i, . . . , 5„\ Sn) 

{Si5j5k)cdidjdk 



y 

6 ^ 

i , J , — ?7i + 1 



W^'^{5o;5i,...,5m;Sm) 



xW^'^{5m;5n 



• 1 ^Ti'i Sjt 



(65) 



Since now the derivatives act only on the second W fac- 
tor, the first W factor factorizes and, after integration over 
. . . , <5m-i, cancels the denominator in the conditional 

probability (44). 

Next, as in MR3, we introduce the notation 



G 



{Sn 



£-£^^{SiS.)S{S,)5{Sk)). 



(66) 



and we expand the correlator as 



{5{Si)5{Sj)5{Sk)} 



E 

),q,r— 



{Sn - Sif (67) 



X {Sn-Sjy{Sn-SkYG[^'''-''\Sn) 



As shown in MR3, in the large mass limit (which is the most 
interesting regime for observing the non-Gaussianities) the 

leading contribution to the halo bias probability is given 
by the term in eq. (67) with p = q = r = 0. We neglect 
subleading contributions, which can be computed with the 
same technique developed in MRS. The discrete sum then 
reduces to {5n)c j k=m+i '^i^jdk- We can now peoceed as 
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in the previous section and perform the shift of variables (19) 
on 5i, for with i = m + 1, ■ ■ ■ ,n — 1. The halo probability 
becomes 

nhaloNGCf^n, Sn\6m, Sm) = 

dSm + 1 ■ ■ ■ dSn-l 

' — OO — oo 

(68) 

where Wmb is the probability density in the space of trajec- 
tories with a moving barrier, so that, in general, 

d5m+l ■ ■ ■ I dSn-lWinb{Sm; Sm + l, ■ ■ ■ ,Sn-l,Sn; Sn) 
-OO J — OO 



t(i) 



t{2) 



(69) 



t(1) 



etc. corrc- 



where, as in eq. (23), the terms H^^, n^,^, ^^^j^, 
spond to the different orders in the expansion of the expo- 
nential in eq. (22), but now for trajectories that start at Sm, 
with Sm. small but finite, rather than at So = 0. 

To compute these expressions we can now use the fol- 
lowing identity, proven in MR3, 



E 

i,j,k=l 
= ^ 

~ dBi 



dSi ... dSn didjdkWmb 



dSn Ilhalo {Sn , Sn \ Sm , Sm) ; 



(70) 



where nhaio(5n, Sn\5m, Sm) is given in eq. (48). The calcula- 
tion of the conditional first-crossing rate J-mh ^G{S„\5m., Sm) 
proceeds by finally taking the derivative with respect to Sn- 



J'mb^GiSn\Sm,Sm) = . ^ e" ^^^T^ 



27r(5„ - S^)3/2 

+ ^—^ — \{B„ - Sm)* 

6V2^{Sn - Sm)"/^ ^ 

— {Bn — Sm)^{Vmn + 2(5'„ — Sm)B'„) 

+2(Bn — Sm)"^ {~{Sn — Sm) + '^mn + {Sn — Sm.)'PmnBn) 
+ {Sn — Sm){Bn — Sm){T'mn -\- 6{Sn — Sm)B'n 

^'PmnBn 2(Sn Sm)'Pmn) {Sn Sm) 
2((S'n Sm)l^mn(l^mn ~\~ {Sn Sm)Bn 

- 4{Sn - Sm)PLn)] e" 

+ — L — , _ " — TTTT [{Bn — Sm)^ — {Bn — Sm^Pmn 



3V27r(S'„ - S'^)B/2 

- {Sn - Sm) + 2Vln] e" , 



(71) 



where the prime stands for derivative with respect to Sn. 
Notice that in the limit of vanishing /nl only the first line of 
the above equation survives. It should also be observed that 
terms proportional to derivatives of the cumulant such as S'i 
(which are subleading compared to the terms proportional 
to S3) can also come from terms with p + q + r ^ 1 in 
eq. (67), that wc have neglected, buth which in principle can 
be computed as in MR3. Expanding again in powers of Sm 
up to Sm , we find the scale-independent NG contribution to 
the Eulerian bias parameters. Normalizing the bispectrum 



Ss{Sn)^^{S'{Sn)), 



(72) 



we find the corrections to the halo bias parameters from NG 

S3{Sn) 



Ab 

+ 
+ 



1 NG 



6S'„(B„ -hP„)2 

Bl(2Sn+Vl+4:SnVnB'n) 
4B„P„(-5'„ +Vn+ SnVnB'n) 
Sn{Sn + Vn{3V„ + SSnB'n - AV^Bn 
S3{Sn)S, 



[3Bt + 2BUrn-2SnB'n) 



3{Bn+Vn)^ 



{{Bn+3Vn){Bn 



WSnV'n)] 
Vn) + Sn) , (73) 



and 
A62NG 



S3{Sn) 



3S'2(S„-h7'„)2 

-h BliVl^ASnVnB'^+d^Sn 

+ Bl{7VnSn - 4Pl - QSlB'n - 



[-iBl + Bi{-2Vn+4.B'nSn) 



ASnV'iB'n) 



+ Bn{4:SnPnB'n - l2SlVuB'n + WSlPnV'n 
- Sn{3Sn + 47?^)) + 2VnSn{Vl - Sn) + 2VlSlB'n 
2S3{Sn) 



S{Bn+Vn) 



{B„{Bn + 3Vn){Bn ~ Vn) " SnVn) , (74) 



where Vn = V{0,Sn). As shown in Fig. 2, these terms cor- 
rect the bias calculated with Gaussian initial conditions by 
a few pcrcents. In the last term wc have expanded for large 
masses and again assumed that the barrier as well as 53 are 
slowly varying. At first one might think that such an ap- 
proximation, altough useful in some cases, would not apply 
to the barrier which corresponds to the the ellipsoidal col- 
lapse, cq. (2). In this case in fact Bst{S„) is given by a con- 
stant plus a term proportional to 5*^ with 7 ~ 0.6 < 1, and 
therefore already its first derivative, which is proportional 
to 5,1 ~^ is large at sufficiently small Sn, and formally even 
diverges as Sn — ^ 0. However one should not forget that, in 
practice, even the largest galaxy clusters than one finds in 
observations, as well as in large-scale A''-body simulations, 
have typical masses smaller than about W^^h~^ Mq which, 
in the standard ACDM cosmology, corresponds to values of 
Sn 5^0.35, see e.g. Fig. 1 of Zentner (2007). Even for such a 
value, which is the smallest we are interested in, the value of 
Bgx(Sn) is just of order 0.3 which means that, in the range 
of masses of interest, the barrier of ellipsoidal collapse can 
be considered as slowly varying. Note that wc have assumed 
a barrier shape which is the same for Gaussian initial con- 
ditions. A justification of this assumption can be found in 
Lam, Sheth & Desjacques (2009). 

For high mass halos, » 1, the scale-independent NG 
contribution to the Eulerian bias parameters are 



ING 



i<S3 [3ozy2-1.6(oi^2)° 



Abo 



S3 



3a' 



3/2 

Sc 



1.6 a" 



(75) 



Notice that the leading term of the first bias coefficient does 
not agree with what found by Smith et al. (2010), who found 
A61NG — —(2/6)53!^^ for high mass halos, even in the limit 
in which wc reduce artificially ourselves to the spherical col- 
lapse case (a — 1) adopted to produce their formula (39). It 
might well be that the discrepancy arises from the fact that 
Smith et al. (2010) deduce the halo bias coefficient from 
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Figure 2. The magnitude of the NG correction (with /nl = 100) 
to the first bias parameter (Eq. (73)), normalized to the halo bias 
without NG (Eq. (51)), as a function of f = 5c/cr. 



the NG halo mass function defined artificially as the Gaus- 
sian Sheth-Tormen mass function multiplied by the ratio 
of the Press-Schecter NG and the Gaussian mass functions. 
On the other hand, had we also considered the variance Sn 
and the skewness .Ss as a function of scale through the long 
wavelength model Sm (k) , we could have obtained as well the 
standard scale-dependent part of the NG bias coefficient, see 
Desjacques et al. (2011). 



4 THE TWO-BARRIER FIRST-CROSSING 
RATE 

Following Lacey & Cole (1993), we now wish to calcu- 
late through the excursion set method the two-barrier first- 
crossing rate. Before launching ourselves into the computa- 
tion of the conditional probability, we pause to make some 
preliminary considerations. In particular, Section 4.2 will 
be devoted to the case of Gaussian initial conditions and 
our aim is there to provide an analytical understanding of 
the problem previously discussed by Bond et al. (1991); the 
case of NG initial conditions will be discussed in Section 
4.3. Bond et al. (1991) showed that the uncorrelated steps 
solution for the two-barrier problem appears to be in better 
agreement with halo conditional mass functions than did 
the solution for correlated steps. Given this, it is natural 
to ask if going from constant to moving barriers alleviates 
the discrepancy. This is excellent motivation for studying 
the two-moving barriers problem. In this case, the uncor- 
related steps problem was solved numerically by Sheth & 
Tormen (2002) who also provided an analytic approximar 
tion, which will will discuss below. Moreover, Sheth & Tor- 
men (2002) also showed that these walks (with uncorrelated 
steps) and their associated forrrmlac for the first crossing dis- 
tribution, were able to provide a rather good description of 
the conditional halo mass function. So the most interesting 
question would be if accounting for correlated steps helps 
in the computation of the two-barrier problem with moving 
barriers. We will not discuss this problem combining Gaus- 
sian initial conditions with non-Markovian evolution of the 
noise for moving barriers. However, we will include the non- 
Markovianity induced by the NG initial conditions. 



4.1 Some preliminaries 

Let us first analyze the problem with only one barrier, _B„. 
An instructive alternative way of understanding why the 
first crossing rate is obtained by taking (minus) the deriva- 
tive with respect to Sn is the following. In order to impose 
that the trajectory makes its first barrier crossing at "time" 
Sn, we require that for all times ^i, . . . Sn-i the trajectory 
stays below Bn, while at Sn it must be above. Therefore the 
quantity that we need is 



/■Bi pB„_i poo 
d6i... / d6n-i / d6nW{6o;Si,...,6 

-oo J — oo J Bji 

We now write 

poo poo pBji 

I dSn — I dSrt — I ddn 

J B-n J ~oo J —OO 



.(76) 



(77) 



and use the fact that 



/Ol /"^n— 1 /"OO 

dSi... / d5n-i / d5nW{5o;5i,...,5 n't Sri) 
-OO J —oo J —OO 

= f dSi... f 

J —oo J —<. 



d5n-lW (So; (5l, . . . , Sn-l\Sn-l) 



dSn-lii-mbi^O', Sn-i; S„-{) . 



The second contribution obtained from eq. (77) is 
/•Si rB„ 

/ d5x... d5nW{5o;5u--.,5n;Sn) 

J —oo J —OO 



(78) 



O — C 



d<5„n^b(<5o;<5„;S„). (79) 
Then we get 

n'-i Sri) 



-f 



n— 1 j Sn—\) 



r 



d(5„nnib(<5o; Sn, Sn) 



= — e-^^ / dSn^mhiSo; Sn; Sn) 

J-oo 
r-B„ 

dSn-i nmb(<5o; <5„-i; Sn-i) + 0{t) , (80) 



where e — {S„ — S„-i). The integral in the last line in the 
continuum limit becomes {B„ — _B„_i)nmb(<5o; Bn; Sn) and, 
since (B„ — B„_i) = 0(e) while nmb(<5o; Sn = Bn; Sn) van- 
ishes as v^, this term is overall 0{e^^^), while the term 
proportional to d/dSn is 0{e). Therefore the transition rate 
per unit time-step e is 



d f"" 

J^mb(S'n) = — ^TT / dSn^inh{So;Sn; Sn) , 

""^^ J-oo 



(81) 



which is the standard result. 

In the two-barrier problem, denoting by B°'{Sni) and 
B''{Sn) the two barriers corresponding to redshift and 
halo mass Mm and Zb and M„, respectively, the relevant 
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quantity is the conditional probability given by the ratio 
between 



/ dSi... dSm-1 i 

J-oa J -oo J 

Bm+1 

dSm+l ■ 

-oo 



dSri 



POO 

X / d5„W{5Q;5i,.. 



(82) 



and 



= / dSl... 

J —oo J —oo 



dSr, 



X / dSm W{5o; 5i,...Sm; Sm) ■ 

■J B% 



(83) 



In eq. (82) the integral over d5m has a lower limit be- 
cause we require that the trajectory crosses above this bar- 
rier. The subsequent evolution can bring it below this barrier 
again, so the integrals over dSm+i, ■ ■ ■ dSn-i have as lower 
limit — oo. Furthermore, the upper integration limit for the 
integrals over dSm, ■ ■ ■ , d6n-i is given by the upper barrier 
B'' because we want to sum only over trajectories that never 
crossed the second barrier 5n at times smaller than Sn. Ob- 
serve that this must be imposed even in the integral over 
dSm- Finally, at S„, the trajectory crosses for the first time 
above Bn, so the corresponding integral runs from B^ to 
-|-oo. The denominator (83) gives the appropriate normaliza- 
tion to the conditional probability, so that in the Markovian 
case, where W factorizes, it cancels against the integrations 
over dSi , . . . dSm in the numerator. Observe also that the 
integral over d6m (both in the numerator and in the denom- 
inator) in the continuum limit can equivalently written as 
in integral from B"^^ and +00, since at time Sm-i the tra- 
jectory was below the lower barrier _B", and the probability 
that in an infinitesimal time step e it jumps above the up- 
per B** vanishes as exp{— (Sj^ — -B^_i)^/(2e)}, so it does 
not contribute to the continuum limit (at least, as long as 
there is a finite separation between the two barriers). 

We now wish to derive a result analogous to eq. (81) 
for the conditional two-barrier first-crossing rate, i.e. for 
the rate at which trajectories first cross the upper barrier 
B'' at S = S„, under the condition that they first crossed 
the lower barrier _B° at 5 = Sm- In the Markovian case 
we can repeat the derivation of eqs. (76)-(81); in this 
case, in fact, W factorizes as W^'^lSo'jSi, . . . ,5n; Sn) = 
W'''^(So;5i,...,5,n;Sm)W'^i5m;5m+i,...,5„;S„ - Sm). 
The fact that the integral over 5m- 1 runs over 5m- 1 ^ Bm-i 
while the integral over 5m runs over 5m S5 Bm implies that 
the factor (27re)"^/^ exp{-(5r„ - 5m-i)^/(2e)} which 
appears in the first W factor (see eq. (27)) becomes, in the 
continuum limit, a Dirac delta which forces Sm to become 
equal to Bm in the second factor W^™, so inside the integral 
we can write 

W^{do;...,S„;S„)= (84) 
W^'^ido; ...,5m; Sm)W^'^{B^;3m+i, . . . , 5„; S„ - Sm) . 

The first W factor, integrated over dSi , - - - , d6m in eq. (82), 
cancels by construction the denominator D, and the integral 



over d6„ in eq. (82) can be treated as in eq. (77). As a result, 
we get again eq. (81), except that So is now replaced by Bm 
and Sn by Sn — Sm. The same argument can be repeated 
in the non-Markovian case. Simply, all derivatives di, when 
acting on the second factor 14^^™ (5m ; 5m+i , ■ . ■ ,Sn', Sn), must 
be evaluated at 5m = Bm (including derivatives di with 
i = m). In this way the dependence on Sm remains only in 
the first W^^™ factor, and the derivation of eqs. (76) (80) 
goes through. 

After all these considerations, we may write the numer- 
ator (82) in general as 



dS'mdS n 



I ' dSi... f 

J —00 J —c 



dSn 



fBl+r PBI 

/ dSm+1 ■■■I dS„ 
J —00 J —00 

[l + f] W^'^{5o;Si,...5m;S'm)\„, _„ 



(85) 



where the differential operator /, which acts on both W 
factors, has the general form 



<=1 i,3=l 
n 

+ c{Si,Sj,Sk)didjdk + .. 



(86) 



i,j,k — l 



and takes into account the effects of the the moving barrier 
and/or the non-Gaussian contributions. The denominator 
(83), as we saw, is given by 



V: 



d 



dSi 



dSn 



(87) 



dSn 

X [l + f]W<''^{So;Si,...Sm;Sm). 

The ratio J\f/I> defining the conditional two-barrier proba- 
bility is therefore proportional to e and the flux J-lnb, which 
is obtained dividing further by e, is therefore equal to the 
ratio of eq. (85) and eq. (87), i.e. 



mb ( Bn , Sn \Bm, Sn 



lim 



(^) 



(88) 



Notice that in the fully Markovian case and constant spheri- 
cal collapse barrier the operator / vanishes. Then factoriza- 
tion applies and one recovers the standard result of Lacey 
& Cole (1993) that the two-barrier conditional probability 
with constant barriers Sb{zb) and 5c(zc) is given by the usual 
first-crossing rate where the initial conditions are such that 
the smoothed density contrast is Sm = Si, at time Sm 



J^5ph(5c, Sn\Sb, 



(5c - 56)e 



-(Sc-Si,)^/(2(Sn-S,n)) 



^{Sn - 5m)3/2 



(89) 



4.2 The two-barrier conditional probability: the 
moving barrier case and Gaussian initial 
conditions 

In the case in which the baxrier threshold is moving, as 
for the Sheth-Tormen barrier, and the initial conditions are 
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Gaussian, the computation of the two-barrier conditional 
probability is very similar to the one we have performed 
in the previous section for the halo bias. The key point is 
to perform, in eqs. (82) and (83), the following shift of the 
integration variables 5i with i ^ m and i ^ n 



Si ^ Si 



b; 



^ pi 



(s, 



1,- 



,"1-1, 



°° nf>,(p) 



^ p\ 

p=i 



(90) 



(where Bm'''' denotes the p-th derivative with respect to S 
of the barrier B'^{S), evaluated in 5* = Sm, and similarly 
B*^^^^ is the p-th derivative of the barrier B*'{S)), while no 
shift is made for the variables Sm and 5„. It is easy to con- 
vice oneself that, exploiting the factorization property (46), 
the contributions to the two-barrier conditional probability 
coming from the random walks starting at 5o = at 5* = 
and crossing the barrier Bm for the first time at Sm cancel 
and one is left with 



b a d /■ " + 1 

W{Bm; Sm+1 ■ ■ ■ ,S„; Sn Sm ) 



dSn 



-(S„-S^)^/{2{S„-Sm)) 



27r(S„ - S„)3/2 



[B'^~B^+P{Sm,Sn)] 



(91) 



This finding coincides with the result obtained by Sheth & 
Tormen (2002); in fact it is valid for any generic moving 
barrier (and therefore also for the ST barrier with a differ- 
ent from unity) and justifies what found in Sheth & Tormen 
(2002) on more rigorous grounds. As was the case for the un- 
conditional mass function, the height of the barrier diverges 
for {Sn — Sm) — >■ oo, so uot all trajectories intersect the sec- 
ond barrier. It seems reasonable (Sheth & Tormen (2002)) 
to associate the fraction of random walks that do not with 
the fraction of the parent halo mass that is not associated 
with bound subclumps. 

Once the two-barrier first-crossing rate is known, we 
may compute the conditional mass function, eq. (3), to study 
how many progenitors at zt are associated with a descendent 
halo mass AIo at Za- In Fig. 3 we show the conditional mass 
functions for the progenitor halos of a descendant halo of 
mass Mo — 10^''h~^MQ at = and look-back times Az = 
{zb — Za) ~ 0.1 and 0.01. We plot the moving barrier result 
Eq. (91) for the ellipsoidal barrier in eq. (2) with a = 0.707, 
and the spherical collapse result Eq. (89) (both with a = 
0.707 and a = 1). 



4.3 The two-barrier conditional probability: the 
moving barrier case and non-Gaussian initial 
conditions 

In the case in which non-Gaussian initial conditions are 
present, we can deal with the problem in the same way we 
have been treating the computation of the halo bias in the 
case of a non-Gaussian theory. Expanding for weak non- 
Gaussianities in the expression (85) brings down the sum 




Figure 3. The conditional mass functions for the progenitor halos 
of a descendant halo of mass Mo = lO^'^h-'^MQ at z = 0, accord- 
ing to Eq. (91) with a = 0.707 (solid red) and to the spherical 
collapse model with a = 0.707 (dashed blue) and with a = 1 (dot- 
ted green). Two different look-back times are shown: Az = 0.1 
(thick lines), and Az = 0.01 (thin hnes). 



i,j,k — l 



{SiSj5k)cdidjdk = ^ {SiSjSk)cdidjdk 

71 

+ ^ {5i5j5k)cdidjdk 

i,j,k — m-'rl 
m n 

+ 3^ ^ {SiSjSk)cdidjdk 

i — l j.fc — m-|-l 

+ 3 ^ ^ {SiSjSk)cdidjdk ■ 



— l fc — m + 1 



(92) 



We now perform the Taylor expansion (67) and retain only 
the leading terms 



{SiSjSk)cdidjdk ~ ^ {Sm)cdidjdk 

— 1 i,j,k — l 

n 

+ ^ {&n)cdidjdk 

z,j,fc— m+1 
m n 

+ 3^ ^ {SmSl)cdidjdk 

i — l _7,fc — m-|-l 

7n 71 

+ 3 ^ ^ {SmS„)cdidjdk ■ 

— l fc— m+1 

(93) 

As Bm < Bn we may further approximate this sum ignoring 
the terms proportional to {Sm)c and (iSm(5,i)c with respect 
to the others 
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{5iSj5k)cdidjdk ~ ^ {S't)cdidjdk 

m n 

+ 3^ ^ (5„.5^)ca,9jafc 



— 1 



i — l j,fe— m + 1 



^ 0.84 
5 0.83 



(94) 



Again, the crucial point is that, when considering a moving 
barrier, the shift of variables (90) does not involve neither 5m 
nor Sn- This allows factorization and cancellations between 
the numerator (85) and the denominator (87). The non- 
Gaussian two-barrier conditional probability J-"mbNG then 
becomes 



J^inhNG{B„, Sn\ 



d 

dSn 



Figure 4. Ratio of the conditional mass functions with and with- 
out NG, for the progenitor halos of a descendant halo of mass 
Mo = IO^S/i-^Mq at ^ = 0, look-back time Az = 0.3, /nl = 50 
and a = 0.707 
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where limb can be read from eq. (69), J^mb(5'm) from eq. (40) 
and 



(95) 3V2^iSn - S'™)5/2 

— {Bn — Bm)'Pmn — {Sn — Sm) + '2Vmn] ■ 



(97) 



ritwo (^n ? Sn j 5m , Sjj 



y^2Tv{S„ — Sm) 



where the prime stands for derivative with respect to Sn- As 
for the third and fourth terms of Eq. (95), they are 



-{S,^-Sn^r/(XS,^-Sm)) 



•^mbNG("BTi5 "^nl-Bm! '5'm) — ^S2.l{Sn ^ Sm) Sn S^, 



1/2 



-{2S° -i„-5m)^/{2{S„-S™)) 
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^ d{B'')^ mb(-B„, Sn\Bm^ Sm) Q-Qa mb('S',r 



2{Bl-5n) 
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(98) 



27r(S„ - S'™)3/2 

2(-B,''i — ^n)^ ^-(2B,'',-g„-f„)2/(2(S„-S„))p2 



(96) 

Notice that repeating with care the steps described in sub- 
section (4.1), in eq. (95) we do not have to differentiate the 
cumulant ((5m5n)c with respect to S'm and Sn- 

Let us set ((5^(5m)c = S2,i{S„, Sm)SnSU^ , where 
52,1 (5*71, Sm) is slowly changing with Sn and Sm (see App. A 
for further details and useful fitting functions). Then, the 
first two terms of (95) read 

•^mbNG("B"' ^^\Bm: Sm) + •^mb'^NG ("^^ ' '^n\Bm^ Sm) 



where is J-mh{Sm) is given by eq. (40) and 

-(IV) 



•^mbNG("B^' '^^i-B"^' '^"^) — ■^<S2.l{Sn: Sm)SnSr^ 



1/2 



■ d(BV)'^dB. 



nib {Bn , Sn \ Bm , Sm ) ■ 



(99) 



Notice that the total NG conditional probability is no 
longer a function of the variables [Sn — Sm) and {Bn — B'^ 
as in the Markovian case. In Fig. 4 we show the ratio of 
the conditional mass functions with and without NG, for 
the progenitor halos of a descendant halo of mass Mo — 
lO^^/i~^M0 at Za = 0, look-back time Az = 0.3, /nl = 50 
and a = 0.707. 
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5 HALO FORMATION TIME PROBABILITY 
WITH NON-GAUSSIANITIES 

In this section we present the results for the probabihty dis- 
tribution of halo formation redshifts with the inclusion of 

non-Gaussianities. 

We follow the convention to define the epoch of forma- 
tion of a halo as the time when the halo contains half of 
its final mass, although the generalization to an arbitrary 
fraction between 1/2 and 1 is straightforward. Let us fix the 
mass Mq and the redshift Za of the descendcnt halo. Then, 
the probability that such halo had a progenitor at redshift 
Zb > Za with mass between Mo/2 and Mo (or, equivalently, 
the probability that the formation redshift is bigger than Zb) 
is given by cq. (3) integrated from Mo/2 to Mq: 

rSh 

P{Zb;Mo,Za) = / dSn 
J So 



Mo 



M{Sn) 

X J-^b(B(«6; Sr,\B{za; So), So) (100) 

where So = S{Mo), Sh = S(Mo/2). From this, it is possible 
to find the probability distribution that the halo of mass Mq 
at Za would have formed between zt, and zt, + dzb'- 



p{zi)dzi = 



dP{zb;Mo,Za 



dzb 



dzb ■ 



(101) 



In the simple case of constant barrier 5c{z) and Gaus- 
sian initial conditions, the distribution of halo formation red- 
shifts is simply given by (see Lacey & Cole (1993)) 



p{zb) = 2a;(0i,)Erfc 



Uj{Zb) 



V2 



doj{zb) 
dzb 



where 



w{zb) 



^S{Mo/2) - S{Mo) 



(102) 



(103) 



if it is assumed a white-noise power spectrum, leading to 
S(M) oc (Lacey & Cole (1993)). 

Now let us introduce the contribution of non- 
Gaussianities. The two barrier conditional probability is 
given by eq. (95). In order to simplify the calculation and 
reach a compact result we make the following assumptions: 
we ignore the cumulant (SmS^), keeping only the leading 
one {6n) (which means we only retain the terms in eq. (97)); 
we consider tSs as a constant and we assume that S{M) 
scales like M~^ , In these approximations, the (normalized) 
probability distribution of formation redshifts in presence of 
non-Gaussianities becomes 



PNG (26) 
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du){zb) 
dzb 



(104) 



This expression is also valid for a spherical collapse model 
with barrier ^/d5c{z) which may be considered as an approx- 
imation to the ellipsoidal model (Giocoli et al. (2007)). In 



order to gain an intuition of the impact of the non-Gaussian 
correction, for Mo = 10^^ MQh~^ and /nl = 50, in the re- 
gions of Zb where the distribution is greater than 0.1, the 
NG term contributes typically less than 10% with respect 
to the Gaussian one. 

In Fig. 5 we show the probability distributions of formar 
tion redshifts of a halo of mass Mo = IQ^^h-^MQ at Za = 0, 
both for the spherical collapse model with constant barrier 
^/aSc and for the ellipsoidal collapse model with barrier Bst 
in eq. (2). These results are obtained by integrating numeri- 
cally eq. (100) using the first-crossing rates in eq. (91) and in 
eqs. (97)- (99) for Gaussian and non-Gaussian initial condi- 
tions, respectively. The inclusion of non-Gaussianities tends 
to shift slightly the distributions towards higher redshifts. 
Furthermore, the results for the spherical collapse with bar- 
rier ^/aSc and ellipsoidal barrier are quite close to each other, 
confirming the suggestions of Giocoli et al. (2007). 

The mass dependence of the variance S{M) enters into 
eq. (100) in an important way. We make use of the numerical 
fit (Al) (solid lines in Fig. 5), but we have verified that us- 
ing a different fit to S{M), like the one of Neistein & Dekel 
(2008), does not change the solid curves appreciably. For 
comparison, we also show the distributions (dashed lines) 
one would obtain by using the simple scaling S{M) oc M~^, 
which is very useful to caxry out analytical calculations and 
it is commonly used in the literature. However, as shown 
in the figure, the use of this simple scaling leads to non- 
negligible differences with respect to a more accurate numer- 
ical fit. Note also that for Zb Za the numerical fit S{M) 
gives values of dp/dzb systematically higher than for the 
case S{M) oc M-\ This point was already noted in Lacey 
& Cole (1993) where it was stressed that using a variance 
with a scaling different from M~^ gives rise to probabilities 
dp/dzb which do not vanish for Zb Za- This is probably 
due to the fact that the particles that the analysis of trajec- 
tories tags as having a given mass M are not really grouped 
into objects of mass M (Lacey & Cole (1993)). 

The analytical prediction (104), which has been found 
for a constant barrier under the assumption S{M) oc M~^ 
and leading NG term, turns out to be in very good agree- 
ment with the numerical result obtained by integrating 
eq. (100) numerically, under the same assumption for S{M). 



6 CONCLUSIONS 

In the excursion set theory the density perturbations depend 

stochastically with the smoothing scale and the computation 
of the halo mass function is mapped into the so-called first- 
passage time problem in the presence of a barrier. Other 
properties of dark matter halos, such as halo bias, accretion 
rate, formation time, merging rate and the formation history 
of halos, can be studied using the excursion set formalism 
by computing the conditional probability with non trivial 
initial conditions and the conditional two-barrier crossing 
rate. 

In this paper we have performed the calculations of such 
conditional probabilities in the presence of a generic mov- 
ing barrier and for both Gaussian and non-Gaussian initial 
conditions. Our generic results can therefore be applied to 
the case of the ellipsoidal collapse where the barrier is mov- 
ing, given by the expression (2), and to the case of the dif- 
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Spherical collapse with barrier '^c 




Ellipsoidal collapse 




Figure 5. Probability distributions for the formation redshift Zf, 
of a halo of mass Mq = 10^^/i~^A/q at Za = 0, in the spherical 
collapse model with barrier y/a&c. (top panel) and in the ellipsoidal 
collapse model (bottom panel). The Gaussian case is shown in 
blue, while the non-Gaussian one, with /nl = 50i is in red. Solid 
lines correspond to the numerical fit to S{M) in eq. (Al) while 
for the dashed lines the simple approximation S{M) oc is 
used. 



fusive barrier discussed in MR2. Our findings include the 
non-Markovianity of the random walks induced by NG. 
Let us summarize the main results of this paper: 

• assuming a sharp filter in momentum space, the first 
two halo bias coefficients for a generic moving barrier 
(eqs. (51) and (52)) and their corrections due to non- 
Gaussianities (eqs. (73) and (74)); 

• the conditional mass function (eq. (3)) for a generic 
moving barrier with Gaussian initial conditions (eq. (91)) 
which reproduces what found in Sheth & Tormen (2002) 
and in presence of non-Gaussianities (eqs. (97), (98), (99)), 
again for a sharp filter in momentum space; 

• the probability distribution of the halo formation time, 
including non-Gaussianities; we have provided numerical re- 
sults for the spherical and the ellipsoidal collapse models 
(Fig. 5) and an analytical approximation (eq. (104)) valid 
for constant barrier, again for a sharp filter in momentum 
space. 

Our calculations have revealed a different scaling of the 
Gaussian halo bias parameter at high halo masses from what 
found in Sheth, Mo & Tormen (2001), our prediction being 
about 20% higher. This seems to go in the right direction 
to fit better the N-body data by Tinker et al. (2010) when 
such a comparison is possible. 



As an application of our findings, it would be interest- 
ing to investigate e.g. the NG halo assembly bias, recently 
discussed in Sheth & Tormen (2004); Gao et al. (2005); 
Reid et al. (2010). 
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APPENDIX A: 
CUMULANTS 



NUMERICAL FITS TO 



We report here a collection of numerical results we have 
found and used throughout the paper. For the window func- 
tion, we have assumed a top-hat in wavenumber space, but 
with a mass-to-smoothing scale relation of a real space top- 
hat filter. The following cosmological parameters from Ko- 
matsu et al. (2010) are used: h = 0.703, Ha = 0.729, = 
0.271, nt = 0.0451, (78 = 0.809. 

The variance S = {5j{) depends on the smoothing scale 
R and, in turn, on the halo mass M (in units of Mc^h"^) 
according to the fitting function 



S{M) 



Co 



C3 
103 



M 



10 

1/6 



2/15 



+ 



C4 



(Al) 



with CO = 7.2 X 10^Cl = 1.2, C2 = -5.8, C3 = 9.5, C4 = 2.6, 
which agrees rather well with the one reported in Neistein 
& Dekel (2008) , who instead make use of a top-hat filter in 
real space. 

For the scale-dependence of 53 = {S^{S))/S'^ we have 
found the following simple fitting formula, computed along 
the lines of Matarrese et al. (2000), with updated cosmo- 
logical parameters: 

53(5') ~ /nL ■ 



SO-'- 

The time- variation of 53 is such that SSySs 
natively, one may define the quantity 



61 (S) 



SS/2 



(A2) 
0.3. Alter- 

(A3) 
- 0.2. This 



which varies more slowly in 5", since Se'i/ei 
quantity is well-fitted by the formula 

ei(S) ~ 2.9 X lO-"* S° Vnl . (A4) 

For the cumulant {5'^{Si)5{S2)), we have found that it 
scales approximately like Si\^S2. Therefore, it is convenient 
to define a slowly-varying <S2,i('S'i, S2) as 



52,1(5*1, 5*2) = 



{SHSi)5{S2)) 



SivS 



(A5) 
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Having fixed 52 = Sir, ~ 0.2, the variance corresponding 
to a halo mass of IO^^Mq/i"^, we have found the fitting 
formula (for S > S15) 



m— 1 Ti — 1 



<: Q . , 2.4x10-^ 



(A6) 



Its very mild dependence on S justifies the assumption made 
in the text, where we consider it a constant. 



APPENDIX B: THE TWO-BARRIER FIRST 
CROSSING RATE WITH A TOP-HAT 
WINDOW FUNCTION IN REAL SPACE 

In this Appendix we would like to extend the computation 
of the two-barrier first crossing rate to the case in which 
the window function for the smoothed density contrast is a 
top-hat in real space. An analogous computation has been 
perfomed by Ma et al. (2010) for the halo bias and we are 
going to use many of the results of the Appendix of that 
paper. The choice of a top-hat filter in real space introduces 
by itself a level of non-Markovianity. The latter is manifest 
in the two-point correlator of the smoothed density contrast 
(MRl) 



{diSi)6iSj)) = min(5i, Sj) + A{S„ S, 



(Bl) 



where A{Si,Sj) = A(Sj,Si) and, for Si ^ Sj, the function 
A{Si, Sj) is well approximated by 

A{Si, Sj) = Aij ~ K ^'(^^^-^') , Si < Sj , (B2) 

Sj 

with K ~ 0.44 (0.35) for a tophat (Gaussian) filter in coor- 
dinate space. The parameter n gives a measure of the non- 
Markovianity of the stochastic process. 

We perform the computation of the two-barrier crossing 
rate assuming a spherical collapse and we will extend them 
to the diffusive barrier model introduced by Maggiore & 
Riotto (2010b) at the end. The barriers are called db and 5c. 

To perform the computation we use the technique dis- 
cussed in detail in MRl. We consider first the numerator in 
eq. (44). The first step is to express the non-Markovian W 
in terms of W^^m^ 

W{5Q;...,5^;Sn) = 

c^W^'^{5o;...,S„;Sn) 
1 " 

+ 2 XI ^ij9idjW^'^{5o; . . . , 5„; S„) . (B3) 

iii = l 

As usual it is convenient to split the sum into various pieces 

m—l 



i,j=i <=1 
_^ n— 1 n— 1 

+- ^ Aijdidj + ^ Aindidn 



i,3 = '^ 



i=m-\-l 



m—l 



(B4) 



i=l j=m-}-l 



j=m-}-l 



where when sums arc from 1 to m the non-Markovian kernel 
has to be thought as a functio of S!^ ■ The goal is to compute 
the numerator (85) and the denominator (87). Consider first 
the contribution from the first line of (B4). Its contribution 
to the numerator in (85) (a part from the time differentation 
with respect to S'^ and Sn) can be written as 



dSi- 



dSrr 



dd-m+i ■ ■ ■ dd„ 



^ m — l m—l 



X W^'^iSo; ...,dm; S'^)W^'^i6m; ...,6n;Sn- Sm) 



_^ m—l m—l 



dSi • • • d6; 

L i,j=l 

xW^'^{So;...,Sm;S'^) 

X / d5m+i---d6nW'^'^{dm;...,d„;Sn-Sm) 

J —00 

dSi--- dSm X ^imdiW'^iSo; ...,Sm; S'm) 

-00 ■_1 



X / d5m+l---d5ndmW^'^{Sm;...,Sn;Sn-Sm). (B5) 

— 00 

The first term is easily dealt with by observing that 



/ 

J — C 



dSm+1 ' ' ' dSnW^"'(Sm; . . . ,Sn;Sn — Sm) 



f 

J — c 



d5n^^^{^m\ (5n; Sn — Sm) • 



(B6) 



Combining this with the contribution coming from the zero- 
th order term W^™^{5o; . . . ,5„; S„) in eq. (B3) and using 
again the factorization property of ly^™, we therefore get 

dS„U^'^ {Sm; S„; S„ - Sm) dSi---d6m 

-00 J —00 

^ m—l m — l 

1 + - ^ Aijdidj + ^ Aimdidn 

i,j=l i=l 

xW''^i6o;...,Sm;S'^) 

/St, 
dSi--- dSm X AimdiW^'^iSo; . . . ,<^m; S^) 
-00 ■_1 



/: 



dSm+l ■ ■ ■ dSndmW'^iSm; ...,Sn;S„- Sm) , (B7) 



where one has to recall that the derivatives dm, acting on 
W^™{5m\ ■ ■ ■ , S„; Sn — Sm) havc to be evaluted at Sm = Sc. 
Notice that the term in brackets gives just the expansion 
to 0{k) of the denominator (87) and therefore it will pro- 
vide the usual Maxkovian tow-baxrier crossing rate (89) . The 
other terms contribute to the numerator (85) as 
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dS'r, 

where 



q2 pSb fSc 

/ d6m / d6n{Na + Ni + + 



,(B8) 



/Ob 
d5i--- d6m-i ^imdiW^'^iSo; ...,5m; S^) 
-OO i=l 

n't Sji Sjri}\ 

J —OO 



d5i---dSm-i / dSm-\-i • ' • d6n-i (BIO) 

-OO — OO 

_^ n — 1 Ti — 1 

- ^ Aijdidj + ^ AindiC 

i,j=m-\-l i=m-\-l 
xW^'^iSo; ...,Sm; S'rr,)W^'^{Sb; ...,Sn;S„- Sm) , 



The term Ni, is given by 

Nb = n<''^{So;5m;S'r,) (B15) 
x[U''^ {5b, Sm; Sn, Sn) + u'''^ {5b, Sm; S„, S„)] , 

where 

d5m+l • • • dSn-l 



Il''^{5b, Sm;5„, S„) = J 
(B9) " X! ^indid„W^{5b;...,5n;Sn-Sm), 



(B16) 



X 

i = m-\-l 

and 



d5i---d5m-i I d5m+i ■ ■ ■ d5„-i 



/d5i---d5m-i I 
-OO J — 

m — 1 



(Bll) 



Nd = I ddi ■ ■ ■ ddm-i I d5m+i ■ ■ ■ d5n-i 



m— 1 n—1 



Aijdidj + ^ Amjdjdn 

_ i=l j=m-\-l j=m-\-l 

xW^'^{5o; ...,5m; S'^)W^{5b; ...,5n;S„- Sm) , 

(B12) 

Notice that in the last term Nd particular attention has to 
be paid on how to reconstruct the derivative with respect to 
S'^ that appears in the numerator (85). Again, we reiterate 
that the second W^'^ has to evaluated at 5m = 5b as well as 
its derivatives with respect to 5m- 

The contribution from Nc vanishes because it contains 
a total derivative 5„ of a quantity that vanishes at 5„ = 5c. 

The term Na is immediately obtained using eqs. (108) 
and (109) of MRl, and is given by 



-V27r 



SbV S'm 



(25b-5mr 



2S' 



+77 



Sm 



Erfc 



25b — 5n 
V^S'm 



X dmli^'^{5m ;5n;Sn Sm)\i^—i^ 

The corresponding flux rate is given by 



(B13) 



J'^Zli^.,Sn\5b,Sm) 



_e 2(S„-Sm) . 

2 {Sn Sm)^^^ 



{Sm -Sr, + {5b - 5cf) e^^Erfc 



' 2*5*77 



(B14) 



n''\5b,S 



(B17) 



;5„,S„)= / d5m+i ■ ■ ■ d5„-i 

J —OO 

n-1 

xi AijdidjW'''^{5b;...,5^;S„- Sm). 

i,j=m-{-l 

The computation of n*"^ and n*^ is quite similar to the com- 
putation of the terms called H""""" and n"'"""-"'"" in MRl, 
and in the continuum limit e — ^ we get 

61 1 f^" 

U''^{db,Sm;dn,Sn) = dnlim - / dSi 

xA{Si, Sn)nr{5b; 5,; Si - Sm)'nr{5c; 5n; 5„ - Si) 

rSn 



= -{5c-5b)d„\{5c-5n) / dSi 



S, 



Sni^Si Sm^'^^'^i^Sn Si)^^'^ 
{Sc — (^c — 



X exp 



2(jSi — S'm) 2(jSn — Si) 



(B18) 



and 



Il^^{5b,Sm;&n,S,^) = Ylm\ I dSi [ dSj 
" J J Si 

xA{S„Sj)nr{Sb;5,;S^- Sm) 
xYinSc; 5c; Sj - S,)Tir{5c; <5„; 5„ - Sj) 
{5c — 5b){5c — 5„) 



-(«c-«(,)V[2(Si-S„)] 



TTV StT 

rSn 



X / dSi 



)3/2- 
e-(«c-<Sn)Vl2(Sn-Sj)] 

5,(S,-Si)V2(S„-5,)3/2 



(B19) 



This can be rewritten as a total derivative with respect to 

5n, as 

Il''^{5b, Sm; 5n, Sn) = 7= {5c — 6b)dn 



7rv27r 



/ dSi 

J Sr. 

Is " '^^^ 



Si 



-{Sc-Sbf/l-2iSi-S„,)] 



{Si - 5m)3/2 

g-(«c-5n)V[2(S„-Sj)] 

5,(5,-S0V2(S„-5,)V2 



(B20) 



The fact that both n''^ and II'''^ can be written as a deriva- 
tive with respect to 5n simplifies considerably the computa- 
tion of the contribution of Nb to the numerator (85) , since 
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we can integrate 9„ by parts, and then we only need to eval- 
uate the integrals in eqs. (B18) and (B20) in 5„ = 5c, which 
can be done analytically, as discussed in MRl. In partic- 
ular the contribution to the numerator from IT*^ vanishes 
because it is a derivative 9„ of a quantity that vanishes at 
6„ = Sc- Following the Appendix of Ma et al. (2010), the 
contribution to the two-barrier first crossing rate from Nb 
can be easily computed to be 



•^sjhi^'^' ^n\Sb, Sm) — ~'gg~ 



sph 



k(5c - 3b) 



(B21) 



The most complicated term is Nd- We get 



Nd = 



277 , e 



— ^ — i± 



- \^2T^bb^S'm^ 

—7rdb{Sb — 5m)Erfc 



2Sb — 5m 









+ — ^^"^{571%', S'm)dm [Sml2 — S^ll] \g 

+ !^dmn^'^{Sm;S'm) [Sml2 - S^Il] |^^^^^ .(B22) 



where 
T2 = 



dS. 



(B23) 



(Sc — Sm.)iSc — Sn) 



' Sm 

X exp 



(Sc — Srn)'^ (Sc — Sn)^ 

2{Sj Sm) 2(5*71 Sj) 



27r 



Sn Sm 



-d„e 



-(2Sc-S,n-S„f /mSn-Sm)] 



(B24) 



J 



dS 



X exp 



{Sc — Sm) {5c — Sn) 

'2{Sj Sm) 2{Sn Sj) 



We are only interested in its value for 5n = Sc as Nd contains 
a total derivative with respect to 5n 



J{5m,5n= 5c) = 



1 g+(«c-«m)V(2S^) 



(B26) 



xErfc 



(<^c — 5m) 



2Sm{Srt Sm) 



The corresponding flux rate is given by 



^iptC-^c, Sr.\5b, Sm) - ^ _ ^^^5/2 



K e 2(s„-s„) 



5b (Sm — Sn + {Sc — Sb)'^) + {Sc — 5b)Sm ( 3 



{5c - db)' 

Sn Sm 



l^Erfc 



y2Sn 



V2Sn {-S'^ + Sm {Sn + {5c - Sb)^)) 



+V^e2^t7sf-I:) J^{Sn - Sm)'^^\5c - 5b) 

V 



xErfc 



Sn Sc — 5b 



+(<5c - 5b) 



Sm y'2{Sn-Sm)_ 

{Sl - 4.SmSl + S^ + SU2Sn + {5c - 5b)^)) 

Sn{Sn Sm) 



^„ J^^^^ (Sn Sm)" 



/2 



^^"^ .V^JA {Sm + {Sc-5b)') 



xErfc 



Sc — 5h 



Sm ^2{Sn - Sn 



(B27) 



The total two-barrier first-crossing rate is finally obtained by 
adding to the usual rate the corrections given by eqs. (B14), 
(B21) and (B27): 



T C u o A _ {5.-Sb)e -i^^-^^> 

:F...{5c,Sn\5b,Sm) - ^(5^ _ Sm)^;2 

+^spl {Sc,Sn 1 5b,Sm)+ J'ipi {5c ,Sn\5b,Sm) 
+J'^^l{Sc,Sn\5b,Sm). 



(B28) 



In the limit Sm <S Sn, the total rate reduces to 

J^sph{5c, Sn\5b) 



1 - K ((5c - Sb)e 2s„ 



2V2^ S^/' 

K 5b 
,3/2 



2nS 



{Sc - 5b) 
S 



27r 5^/' 

{Sc-Sb)^ ^ 

2Sn 
2 



(B29) 



which reproduces the result of Ma et al. (2010). In the case 
of a moving barrier, we expect that it is a good approxima^ 
tion (Giocoli et al. (2007)) to simply replace the constant 
barrier Sc with ^/aSc, and k with an (Maggiore & Riotto 
(2010b); Ma et al. (2010)). If so, we obtain 



(B25) J^mh{Sc,Sn\5b) = 



1 - OK ^/a{5c - 5b)e 2i,'„ 

^3/2 



2n 



, 3/2 K {5c - 5b) p a{Sc - 5b)'^ 
2^/2^ sT V 25„ 



,3/2 



K Sb 



2^ s'i^^ 



a{5c - 5b)^ 



Sn 



e 2s;r- 



(B30) 



Instead, the case of the diffusive barrier is obtained from 
eq. (B29) by simply sending Sn into S„/a and k into an, 
that is 



J^d\f{5c,Sn\5b) 



_a(£c-£br 

l-g/t 3/2 {5c - Sb)e 



+^^izi!ir(o 



2^2^ 53/2 

a^/^K 5b 
2^ 5^/2 



27r sT 
a{5c — 5b)'^ 

2Sn 



a{Sc - Sb)'^ 



e 2s:f- 



(B31) 



We expect these results to be correct in the high mass limit. 
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